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ABSTRACT 

The mass and structural evolution of massive galaxies is one of the hottest topics in galaxy formation. This 
is because it may reveal invaluable insights into the still debated evolutionary processes governing the growth 
and assembly of spheroids. However, direct comparison between models and observations is usually prevented 
by the so-called §progenitor biasT, i.e., new galaxies entering the observational selection at later epochs, thus 
eluding a precise study of how pre-existing galaxies actually evolve in size. To limit this effect, we here gather 
data on high-redshift brightest group and cluster galaxies, evolve their (mean) host halo masses down to z = 0 
along their main progenitors, and assign as their “descendants” local SDSS central galaxies matched in host 
halo mass. At face value, the comparison between high redshift and local data suggests a noticeable increase in 
stellar mass of a factor of > 2 since z ~ 1, and of > 2.5 in mean effective radius. We then compare the inferred 
stellar mass and size growth with those predicted by hierarchical models for central galaxies, selected at high 
redshifts to closely match the halo and stellar mass bins as in the data. Only hierarchical models characterized 
by very limited satellite stellar stripping and parabolic orbits are capable of broadly reproducing the stellar mass 
and size increase of a factor ^2-4 observed in cluster galaxies since z 1. The predicted, average (major) 
merger rate since z ~ 1 is in good agreement with the latest observational estimates. 

Subject headings: cosmology: theory - galaxies: statistics - galaxies: evolution 


1. INTRODUCTION 

The size evolution of massive, bulge-dominated galaxies 
has now tur ned into one of the hottest topics in galaxy for¬ 
mation te.g.. lTruiillo et al.]l2007h . While there are promising 
ideas on how this evolution actually occurred, its exact nature 
is still debated. 

One of the main limitations is that samples usually study 
galaxy evolution at fixed stellar mass, a procedure that in¬ 
evitably includes contributions from pre-existing galaxies in 
the sample evolving in time, as well as new galaxies en¬ 
tering the sel ection at later times . As emphasized by sev¬ 
eral authors (iHopkins et al.l l2009t ICarollo et al.l 120131) . the 
impact by newcomers on the global structural evolution of 
massive galaxies, usually termed as “progeni tor bias” (e.g., 
Ivan Dokkum & Franxlll996l : ISaglia et alJlioiOh . might in fact 

* School of Physics and Astronomy, University of Southampton, 
Southampton S017 IBJ, UK; F.Shankar@soton.ac.uk 

^ Jet Propulsion Laboratory, California Institute of Technology, MS 169- 
234, Pasadena, CA 91109, USA 

^ Department of Astronomy, California Institute of Technology, MS 249- 
17, Pasadena, CA 91125, USA 

Centre for Astrophysics, Cosmology & Gravitation, Department of 
Mathematics & Applied Mathematics, University of Cape Town, Cape Town 
7701, South Africa 

^ Department of Physics and Astronomy, University of Victoria, Victoria, 
British Columbia, V8P lAl, Canada 

® CITA National Fellow 

^ Department of Physics and Astronomy, California State Polytechnic Uni¬ 
versity, Pomona, CA 91768, USA 

* GEPI, Observatoire de Paris, CNRS, Universite Paris Diderot, Paris Sci¬ 
ences et Lettres (PSL), 61, Avenue de TObservatoire 75014, Paris France 

® Department of Physics and Astronomy, University of Pennsylvania, 209 
South 33rd St, Philadelphia, PA 19104 

Universite Paris Denis Diderot, Universite Paris Sorbonne Cite, 75205 
Paris Cedex 13, France 

” California Institute of Technology, Pasadena, CA 91125, USA 
International Center for Theoretical Physics, 34151 Trieste, Italy 


even be the dominant factor controlling the inferred structural 
evolution. There is indeed evidence for a significant increase 
in the number density of massive ear ly-type galaxies with 
cosmic time at least since z 5, 2 (e.g., iBuitrago et al.l l2013l 
iHuertas-Companv et al.ll2013M . 

To overcome progenitor bias, there have been mainly 
three methods proposed in the literature. By assuming 
number density conserva tion in the stellar mass function, 
Ivan Dokkum et all (120101) attempted, on statistical grounds, 
to trace back what is the expected average evolution along 
the progenitors. This methodology , now widely applied 
from intermediate-mass (iPatel et al.l 120131: iPanovich et alJ 
20141), to very massive galaxies (e.g., lMarchesini et al. 2014t 
Ownsworth et al.ir2014t) . clearly presents limitations, such as 
the mass-dependent specific star formation rate distribution 
and the impact of mergers, which could potentially alter the 
galaxy rank-ordering in time (e.g.. Ivan Dokkum et aLlUOlOl : 
iBehroozi et al.l 1201^ . Another p ossibility is based o n ap¬ 
plying some “age filtering” (e.g.. ICimatti et all l2012l) . i.e., 
by considering only galaxies having ages compatible with 
the cosmic time pa ssed from high to low redshifts (see also 
iKeating et al.ll201^ . Both of the former techniques are ac¬ 
ceptable ones to grasp some initial insights into the admit¬ 
tedly challenging issue of empirically defining patterns along 
the elusive progenitor-de scendant evolutionary tracks. _ 

A third approach (e.g.. lLidman et al.ll2012llLin et al.ll2013h . 
which is the one followed, refined, and extended in this work, 
has been to focus only on a specific class of galaxies, namely 
the Brightest Cluster Galaxies (BCGs). The idea is to start at 
high redshift with BCGs matched in stellar and halo mass to a 
given observational sample, calculate their host cluster mass 
growth down to z = 0 as expected from hierarchical models of 
structure formation, and finally infer their average stellar mass 
growth by selecting as descendants those galaxies matched in 
present-day cluster mass via some empirical stellar mass-halo 
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mass relation. 

The aim of this work is to extend this methodology to 
specifically include sizes. To this purpose, we perform a close 
comparison between state-of-the-art hierarchical models and 
a number of measurements of the mass and size of central 
galaxies in low and high-redshift groups and clusters. 

This is one of the few attempts to probe mass and structural 
evolution avoiding progenitor bias in both the models and the 
data. Moreover, BC Gs are the ideal systems f or probing struc¬ 
tural evolution (e.g. jEdwards & Pattonl20f^ . Since the most 
massive dark matter halos assembled their mass recently by 
major mergers, it is likely that the structural parameters of the 
most massive galaxies at their centers are also affected. 

2. DATA 

This work make s use of a cluster sample drawn by 
iLidman et al.l (1201 2b in the redshift interval 0.85 < z < 1.63 
from the SpARCS survey (iMuzzin et al.ll200^ IWilson et alJ 
l2009t iMuzzin et al.l 1201 21) . All 12 clusters were discovered 
searching for overdensities of red galaxies using a combi¬ 
nation of Spitzer IRAC and ground-based z-band images. 
The BCG sizes adopted here are derived by Rettura et al. 
(2015) using GIM2D, a fitting algorithm for parameterized 
two-dimension al modeling of surface brightness distribution 
(ISimardlll998h. to fit each galaxy light distribution by adopt- 
ing a lSersicI d 1963b profile. Halo masses in Sp ARCS, inferred 
from line-of-sight veloci ty dispersions (e.g., iDemarco et alJ 
l2010t lAscaso et al.ll2014b . span two decades in cluster mass 
0.043 X 10^^ < Mhaio/4^o < 2.629 x lO'^ (Wilson et al., in 
prep.). 

We complement th is sample with data on groups at 
z ^ 1 from COSMOS dF i noguenov et al.ll2007b iGeorge et al.l 
1201 U[Huertas-Companv et al.l2013bl and references therein). 
These have halo mass measurements derived from weak 
lensing, and galaxies by single Sersic profiles. The galax¬ 
ies used here have been strictly selected to have a prob¬ 
ability of being massive ellipticals with PROBEll> 0.5, 
based on the Bayesian automa ted morphological classifier by 
[Huertas-Companv et al.l (1201 lb . though we stress that none of 
our results re lies on such a specific se l ection. As extensively 
discussed by Huertas-Companv et al.l (1201 3bb . the resulting 
subsample of galaxies ha ve high Sersic indices with n > 5. 

At z < 0.2 we use the iBernardi et alJ (12013112014b galaxy 
sample extracted from the lMeert et alJ (1201 5b database which 
has new and updated photometric measurements of the Sloan 
Digital Sky Survey (SDSS)-DR7 main galaxy sample (specif¬ 
ically improved sky subtraction and different profile fitting). 
The host halos mass es have been assigne d by match i ng the 
IBernardi et al.l (12013b galaxy sample to the I Yang et al.l (l2007b 
halo catalog. In particular, we recomputed the total stellar 
mass of each group (M* ) in the Yang et al. catalog (using the 
Meert et al. updated photometry), and rank-ordered the halo 
mass to this updated M*^. 

For completeness, we also assign galaxies to halo masses 
using the recent dire ct abundance matchin g results from 
IShankar et al.l (1201 4bb . based on the revised lBernardi et al.l 
(1201 3 I) stellar mass function and a n egligible scatter in stel¬ 
lar rn ass at fixed halo mass (see also iRodrfguez-Puebla et all 
I 2 OI 5 L and references therein). 


• Within our numerical dark matter catalogs, we select 
those halos that fall in the same bin of halo mass of our 
z 1 data (Sect.|2]l. 

• Among these halos we select those that host a central 
galaxy within the range of stellar masses measured in 
the data. 

• We follow the merger trees of the latter subsample of 
halos along their main progenitor branches down to z 

0 . 

• In the local Universe we then select within the SDSS 
catalog (cfr. Sect. |2]), those central galaxies that share 
the same mean host halo mass as the descendants of our 
z 1 halos. 

• The comparison between “descendent” and “progeni¬ 
tor” systems linked by their main progenitors branches, 
informs us on the possible mean growth histories un¬ 
dergone by our z 1 central galaxies. 

• We finally compare the predictions of state-of-the-art 
semi-empirical and semi-analytic hierarchical galaxy 
evolution models for the growth in mean stellar mass 
and effective radius with the ones inferred by our com¬ 
bined z 1 and z < 0.3 data. 


4. MODELS 


In the following, we take as a reference the s emi-empirical 
model d eveloped and extensively discussed byfShank^^^T 
( 2014ab . b ased on the general lines exposed bv lHopkins et al. 
(l2009h and IZavala et al.l (l2012b . We here recall its main ingre¬ 
dients. 

The model is built on top of numerical merger trees ex - 
tracted from the Millennium I Simulation'^ (ISr)ringelll2005b . 
Central galaxies, associated with the main progenitor branch 
of the tree, are initially defined as disks at early epochs. All 
their mass and structural properties are assigned via empir¬ 
ical, time-dependent relations, mor e significantly the st ellar 
mass-halo mass relation, taken from iMoster et al.l (1201 3b . the 
disk radius-stellar mass relation, which we write as (e.g. , 

IShen et ai1l2003l : iHopkins et alJl200^ ISomerville et al.ll2008h 


^disk ~ 


Ro 


(1+z) 


-M 

A 4 Ki ct 


-0.4^’^^star 


id- 




3.98 X IO'OM© 


p-k 


( 1 ) 


with Ro=0.1, k = 0.14, p = 0.39, the specifi c star formation 
rate-stellar mass relation (iKarim et al.ll201 lb , and the mass- 
and ti me-dependent gas fraction as given by IStewart et alJ 
( 1200 ^ . 

Central galaxies are re-initialized at each timestep during 
the evolution and can gradually transform their morphology 
via mergers with satellites. Satellite galaxies are those asso¬ 
ciated with each dark matter branch merging with the main 
progenitor. When galaxies become satellites in larger haloes, 
they are assigned a dynamical friction time-s cale for final co¬ 
alesce nce with the central galaxy taken from lMcCavana et alJ 
(l2012b . They are also assigned all the mass and structural 
properties of a central galaxy in a typical halo, within the 
Monte Carlo catalog, of the same mass at the time of infall. 


3. METHODOLOGY 

Our methodology can be summarized in the following 
steps: 


By adapting the lStiankar et al.l 12014ali basic routines to also ran on an¬ 
alytic merger trees, mainly characterized by lower values of trg = 0.8, instead 
of erg = 0.9 of the Millennium Simulation, we have checked that our main 
results do not depend on the exact underlying cosmology. 
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Fig. 1.— Left'. Predicted mass evolution in the total, bulge, disk, and gas components, as labeled, for a central galaxy in a host dark matter halo of mass 
Mhaio ~ 1.5 X IO'^Mq at z = 0. Right: Comesponding evolution in the half-mass radius of the bulge and disk components; the solid, black line marks the 
mass-average of the full . The galaxy undergoes a major merger at z ~ 0.9 that converts it from a gas-rich disk to a gas-poor elliptical. 


The evolution of satellites after infall takes into account extra 
mass and size growth in the disks from the residual gas, and 
may also include gas and stellar stripping in some models. 

When a merger between two progenitor galaxies of masses 
Ml and M2 is major (M2/Mi > 0.3), it is assumed that a single 
elliptical galaxy is produced, and any gas present in the disks 
of the merging galaxies is converted into stars in a burst. In 
a minor merger (M 2 /Mi < 0.3), all the stars of the accreted 
satellite are added to the bulge of the central galaxy, while 
the gas is added to the main gas disk. New stars are also 
fo rmed in a merger follow ing the collisional starburst model 
bv ISomerville et aO (1200 Ih . where a fraction 

of the cold gas of the progenitors is converted into stars. 

The size of the remnant is computed from the energy con¬ 
servation between the sum of the self-bindi ng energies of th e 
progenitor galaxies, and that of the remnant (ICole et al.l2000h . 

(M1+M2)" _M\ ^ Ml ^ /o,b M1M2 
Rnevi R\ Rl C R1+R2' 
where M,-, are, respectively, the total masses and half-mass 
radii of the merging galaxies, and the form factor c ~ 0.5. 
The factor /orb parameterizes the average orbital energy of the 
merging systems, ranging from zero to unity for parabolic and 
circular orbits, respectively. Eq. |3] can further be adapted to 


include any type of 

energy loss such as gas dissipation (e.g.. 

Hopkins et al.l2009l; 

Coyington et al.l201 lllPorter et al.l2014l; 

Shankar et al.ll2014a 

). 


Our semi-empirical model does not include other “in- 
situ” processes, such as major disk instabilities (e.g., 
iBower et al.l l200 6|l or disk regrowth after a major merger 
(e.g., iHammer et al.l l2009l) . We note, however, that ex¬ 
tensive semi-analytic models have proven disk regrowth to 
play a significant role mainly in intermed iate-mass galaxies 
(iDe Lucia et al.ll2()TTt IWilman et al.ll201^ . below the stellar 
mass threshold considered in this work. We also note that our 
massive ellipticals, remnants of a major merger, are usually 
devoid of any left-over gas, which is mostly consumed in the 
star formation burst (Eq.lU see also Eig.[TJ. In some sporadic 
cases, some later gas-rich minor mergers may be able to re¬ 
grow a light stellar disk, but this has negligible impact on any 
of our results. 

Fig- El shows an example of the mass and structural evolu¬ 
tion of a central galaxy in one of the most massive halos in our 


catalog, with mass Mhaio ~ 1 .2 x 10 '^Mq at z = 0. The galaxy 
starts as a disk at z ^ 3 (left panel), with a small central bulge 
built via minor mergers, reduces its gas reservoir, and finally 
fully transforms its morphology, turning, after a major merger 
at z ~ 1.5, into an elliptical with no residual gas. The sub¬ 
sequent evolution is mainly driven by dry minor and major 
mergers, with its half-mass radius increasing by a factor of 
~ 10 since its first major-merger event (right panel). 

For completeness, in the following we will also compare 
our data with the publ icly available‘^h i erarch ical models by 
IGuo et af] (1201 Ih and iHenriaues et al.l (l2014l) . Specifically, 
we have analyzed several thousands of merger trees extracted 
from their online catalogs, i.e., those contributing to all central 
galaxies with final stellar mass Mstar ^ 3 x 10"Mq. We dis¬ 
cuss the rele vant specifications and predictions of such models 
in Sect. 15.31 

5. RESULTS 

5.1. The overall size evolution ofBCGs 

We start in Fig. 3bv re porting the observational results by 
Ivan Dokkum et ah 720101 cyan region). They fitted the mean 
size evolution of z = 0 galaxies with Mstar ^ 2 x lO" Mq along 
the “progenitors”, identified via number density conservation 
(Sect [T]i, as Re(z) oc (1 -Fz)''^ (see also, e.g., ICimatti et alJ 
l2012h . Their results are compared to the main s emi-empirical 
model realizations, as labeled, developed by IShankar et alJ 
(l2014ah for galaxies selected to have the same cut in stellar 
mass at z = 0, and followed in their mean half-mass radius 
growth along their main progenitor branches. 

All models point to an increase in median size of a factor of 
< 2.5 since z 1, in broad agreement with the data. However, 
not all models perform equally well. Models characterized by 
high orbital energies (/orb ^ 0.5; dashed, orange lines) have 
already been sho wn to predict shallower size -mass relations 
compared to data (IShankar et al.ll201^l2014ah . Fig.|2]extends 
these results by showing that the latter type of models also 
induce a weaker size evolution with respec t to the data. As 
already pointed out by IShankar et alJ (1201 3l) . high-resolution 
cosmological simulations aimed at studying the orbital pa¬ 
rameters of major mergers of cold dark matter halos, find in- 
deed that a large fraction o f encounters are nearly parabolic 
(iKhochfar & Burkerlll2006h . 

Stellar stripping in the infalling satellites is included in one 
variation of the reference model by following the empirically 

http://www.g-vo.org/MyMillennium3 
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Fig. 2.— Predicted mean size evolution, normalized to the value at z = 0 for different models, as labelled, for galaxies that today have a stellai' mass 
logMstar > 11.3 Mq. The left panel considers averages all galaxies in all environments, while the right panel only considers galaxies in cluster environments, 
with logMhaio > 14.5M0, which have a stronger size evolution. Models with significant str ipping in the infalling sat ellites, larger orbital energies, and lower 
dynamical frictions timescales are disfavoured by the derived obseiwational estimates by, e.g.. lvan Dokkum et alj flOig) {cyan region). 


motivated mass-dependent stell ar and gas stripping, parame¬ 
terized as jCattaneo et al.ll2Ql lb 

Fstrip = (l-r?)", (4) 

with T = fdf/fdyn the ratio between the dynamic al friction and 
dynami cal timescales. It was discussed by IShankar et alJ 
(l2014ah that a yalue of t = 0.25 allows to faithfully match 
the local iMoster et al.l (12013h stellar mass-halo mass relation. 
Without stripping, the resulting z = 0 stellar mass-halo mass 
relation would end-up being steeper, in better agreement with 
other, more recent deter minations of th e stella r mass-halo 
mass relation based on the lBernardi et al.l (12013h stellar mass 
funct ion and BCG measureme nts (see, e g.. iKraytsoy et all 
1201^ andlShankar et al.ll2014bll . 

In ishankar et al.l (l2014al) we showed that models with stel¬ 
lar stripping tend to fall short in matching the local size dis¬ 
tributions especially at the high-mass end, mainly due to the 
merging satellites turning too compact to promote efficient 
size growth. Here we extend such result to higher redshifts, 
emphasizing that the same models with stellar stripping (triple 
dot-dashed, purple lines in Fig.|2]) tend, as expected, to gener¬ 
ate a size evolution too weak with respect to what empirically 
inferred. 

Models with significantly lower dynamical friction 
timescales, specifically 1 /3 of the ones used in the reference 
model (dotted, cyan lines), increase the number of mergers 
on the central galaxy and have been shown to be disfavored 
with respect to the local data in terms of their large intrinsic 
scatter and more pr onounced environmental dependence 
dShankar et al.ll20l'4ih . These models also tend to produce a 
faster growth with respect to observational results, especially 
at z > 1, though per se deep data are still too sparse to be able 
to firmly rule out this model, as further discussed below. Gas 
dissipation in major mergers (long dashed, red lines) does 
not signi ficantly alter the size evolution along the progenitors 
(see also lShankar et al.|l2013l) . 

In the following, we will adopt the /orb = 0, no-stripping, 
and no-dissipation semi-empirical mod el as our reference one , 
as it performs better against local data (IShankar et al.ll2014^ . 
and predicts reasonable degrees of size evolution. 

Probing size evolution at fixed stellar mass, besides the ef¬ 
fects of progenitor b ias, could also mix different environments 
at dif ferent epochs (IValentinuzzi et al] 120101 : [Poggianti et alJ 
12013l) . preventing a clear understanding of the true role of en¬ 


vironment in size evolution. In fact, as evident from the right 
panel of Fig. |2l the subsample of central galaxies specifically 
hosted by massive clusters above Mhaio ^ 6 x at z = 0, 

share quite similar mean evolutionary histories with respect 
to all galaxies of similar stellar mass, at least at z 1, while 
they tend to have an accelerated evolution at z > 1. Although 
environment may not play a dominant role in determining 
the degree of size ev olution (iHuertas-Comnanv et al.ll^l3al : 
IShankar et al.l l20l'43) . it may still accelerate the growth of 
galaxies living within the depth of l arge potential wells 
ypassbender et alJ[2014t iMei et all 12014t) . Signs of an accel¬ 
erated evolution in (proto)clusters is now becoming more and 
more evident in high redshift data fro m different groups (e.g., 
iLani et all 1201 31: iDelave et alj|20l4 and references therein), 
and we will reserve separate work on this and other related 
issues. 

5.2. The mean and median growth of host halos 

The comparison between models and data attempted in the 
previous section is only a very approximate one. The tech¬ 
nique of number conservation is in fact only a rough approxi¬ 
mation for tr acking progenitor . Indeed, we verified with the 
aid of the full lGuo et al.l(l201 lb merger tree catalog of massive 
galaxies, that selecting at fixed number density tends to induce 
a somewhat weaker size evolution than the true one. We thus 
proceed below with a more accurate model-to-data compari¬ 
son based on exploiting the high and low-z BCG data listed in 
Sect.|2] following the methodology described in Sect. [3] 

However, before proceeding with such a comparison be¬ 
tween theory and observations, it is essential to point out the 
implied differences between median and mean dark matter 
halo growths. Fig. [^reports the mean (orange, long-dashed) 
and median (purple, long-dashed) halo mass growth of halos 
with mass logMhaio/A^© = 13.7 (left) and logMhaio/440 = 14.5 
(right) from z = 1.2 down to z = 0. These curves are obtained 
by integratin g the respective mean and median accretion rates 
presented bv iFakhouri et al.l (l2010h . derived from direct fits to 
the Millennium simulations. Mean accretion rates are at all 
redshifts a factor of about two higher than median ones, thus 
implying a much more prominent mass growth. 

The solid, red lines in Fig. [3mark instead the evolution of 
the (minimum) halo mass for which the halo cumulative num¬ 
ber density at all subsequent times equals the cumulative num¬ 
ber density above a halo mass of logMhaio/MQ > 13.7,14.5 
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Fig. 3. — Analytic fits bv IFakhouri et 3] {MS) to the mean (orange, long-dashed) and median (purple, long-dashed) accretion rates for halos with mass of 
logMi,aio/A/Q = 13.7 (left) and logMi,aio/MQ = 14.5 (right) at z = 1.2. The solid, red lines show the evolution of the (minimum) halo mass for which the number 
density equals the one at z = 1.2 of halos with mass logM[,aio/4f0 = 13.7,14.5. Mean accretion rates more closely track the mass evolution corresponding to the 
conservation of cumulative number density. 


at z = 1.2 (left and right panel, respectively). The latter evolu¬ 
tionary tracks nicely follow the mean accretion rates down to 
z ~ 0, thus proving that mean accretion rates are better suited 
to track the average mass evolution of a subsample of halos 
from high to low redshifts. In the following we will thus al¬ 
ways adopt mean quantities to describe the evolution of galac¬ 
tic properties. 

5.3. Specific comparisons between model predictions and 
observational data 

Fig. m compares the halo mass, stellar mass, and size evo¬ 
lutionary tracks predicted by hierarchical models to observa¬ 
tional data. The filled circles refer to high-redshift data from 
COSMOS at z ~ 0.9 (upper panels) and SpARCS at z ~ 1.2 
(lower panels). In passing, we note that the z ~ 1 stellar and 
host halo masses in our refere nce data are, within the nomi¬ 
nal scatter of ~ 0.15 dex (e.g.. lRodrfguez-Puebla et aTll2015l: 
IShankar et al.llMl4bh . in line with what predic ted by abun¬ 
dance matching models at these redshifts le.g.. iMoster et al.l 
l2010HShankar et al.ll2014bh . 

Local data from SDSS are marked by open symbols. The 
stellar masses (middle panels) and corresponding effe ctive 
radii (right panels) in ferred from the lYang et aP (l2007h halo 
catalog matched to the lBernardi et al.l (1201^ SPSS DR7 sam¬ 
ple, are shown with open s quares, while those de rived from 
direct abundance matching (IShankar et al.ll2014bl) are shown 
with open triangles. The latter, in line w ith other groups 
(iGonzalez et al.l 20131: iKravtsov et al.ll20r^ , predict at fixed 
halo mass stellar masses (and thus rel ated sizes) in very good 
agreement with those derived from the lYang et al.l (l2007h halo 
catalog for halo masses Mhaio ^ 2 x IO'^^Mq (upper panels). 
For mor e massive host halos , revised abundance m atching 
models (IKravtsov et al.l 12014t IShankar et al.l 12014bl) predict 
stellar masses systematically higher by a factor of ^ 2, and 
consequently effective radii higher by a factor of ~ 1.6. 

At face value, the simple comparison between the high-z 
(filled circles) and local SDSS data (open symbols), irrespec¬ 
tive of the exact stellar mass-halo mass mapping, suggests a 
significant growth of a factor of > 2 in mean stellar mass since 
z ^ 1.0, and a comparable increase of a factor of > 2.5 in 
mean effective radius for central galaxies in both the group 
and cluster environments. These trends support an average 
growth rate of cx with a > 1.3, thus stronger than 

a model based on only major mergers (a = 1), indicative of 
the possibly significant role played by minor mergers in the 


global structural growth. 

We now proceed in checking if the latter, observationally- 
based findings reported in Fig. |4] are in line with the predic¬ 
tions of state-of-the-art, galaxy evolution hierarchical models. 
Fig. |4]is one of the first examples in the literature in which a 
careful comparison is carried out between models and data 
to probe the parallel evolution of central, massive galaxies in 
halo mass (left panels), stellar mass (central panels), and ef¬ 
fective radius (right panels). In this section we detail the out¬ 
puts of the semi-empirical model, while we devote the next 
section to the outcomes of advanced semi-analytic models. As 
outlined in Sect.[3l for a careful comparison between models 
and data we first select in our mock catalogs the progenitor 
host halos that at z ~ 0. 9 and z ~ 1.2 have masses compar able 
to the clusters used by [Huertas-Comnanv et al.l (1201 3bl top 
panels) and Rettura et al. (2015; bottom panels), respectively. 
Among these progenitors we then select those halos that have 
a central galaxy with stellar mass comparable to those inferred 
for the central group and cluster galaxies in the reference data. 
The main progenitor branches competing to the selected ha¬ 
los are then followed forward in time, and at each timestep we 
compute their mean host masses, stellar masses, and effective 
radii. The theoretical mean and 1 a uncertainties are shown 
as solid lines and gray areas, respectively. 

We first note that the mean halo masses of our selected ha¬ 
los (so lid lines in the left pan els of Fig.|4]i evolve as predicted 
by the lFakhouri et al.| (1201 Oh analytic fits, within an accuracy 
of < 0.05-0.1 dex (dot-dashed, purple lines, left panels). The 
mean stellar mass of their central galaxies instead is predicted 
in the reference semi-empirical model (solid lines in the left 
panels of Fig. S to evolve by a factor of ~ 2 since z < 1, 
irrespective of the exact environment, either a group or a clus¬ 
ter. This growth is broadly consistent with both observational 
data sets at z ~ 1 and z 0.1, especially for lower-mass halos 
(upper panels), where the local constraints from abundance 
matching and group catalogs are in very close agreement. 

The size growth is analogously similar in the two en¬ 
vironments, possibly predicted to be marginally weaker in 
the group environment. To appropriately co mpare with the 
data we follow the same approach pursued in IShankar et al.l 
(l2014al) . We first assign to each mock galaxy a Sersic in- 
de x n from the empirical relation calibrated at z = 

0 (iHuertas-Comnanv et al.l 1^13bl) . and then computed the 
correspon ding 3D-to-2D size correc tion from the tabulated 
values of iPrugniel & SimienI (Il997l) . We assume constant 
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Fig. 4. — Mean halo mass (left), stellar mass (middle), and projected half-mass/light radius (right) evolutions predicted by the reference semi-empirical model 
(s olid lines with gra y bands), the semi-empirical model with lower dynamical friction timescales (triple dot-das hed, orange li nes), and the semi-analytic model 
of IGuo et alj (red long-dashed and dotted lines). The upper panels compare with the data by COSMOS jHuertas-Companv et al.ll2()13bl) (filled circles) 

and SDSS local data, while the lower ones with the recent d ata by SpARCS (R ettura et al. 2015, in prep.) a nd SPSS. Local SPSS data on stellar masses and sizes 
(middle and right panels, respectively) extracted from the dYang et ai] |2007|) halo catalog matc hed to the IBemardi et al.l d2014l) SDSS DR7 sample are shown 
with open squares, while those derived from direct abundance matching Ishankar et al.ll2014bll are shown with open triangles. Models broadly reproduce the 
data, especially in group-sized halos. Lower dynamical friction timescale models tend to provide faster growth tracks, in apparent better agreement with massive 
cluster environments. 


Sersic distribution for bulge-dominated galaxies sin ce z ~ 1 
dHuertas-Companv et al.ll^l3bl: lAscaso et al.ll2014t) . though 
we verihed tha t even allowing for signih cant evolution, e.g., 
n oc (1 -Fz)“*^ ®^ (Ivan Dokkum et al]| 201 0 h. yields quite similar 
results, with only a slightly faster size evolution at z > 0.5 in 
the latter case, but still well consistent within the 1 a uncer¬ 
tainties. Our model result s are in broad agreem ent with the 
numerical outputs by, e.g., iLaporte et al.l (1201 3l) . specihcally 
targeting the evo lution of BCGs, and also with those by the 
lOser et all (l2012h . though the latter were mainly focused on 
lower mass halos than the ones considered here. 

We note that at the massive cluster scale, for host halos 
above Mhaio ^ 5 x 1 O'"^M 0 , the reference model tends to fall 
short in fully reproducing the local mean stellar mass and 
size as predicted by the most recent abundance matching re¬ 
sults (open triangles). A semi-empirical model characterized 
by lower dynamical friction timescales (triple dot-dashed, or¬ 
ange lines), tends to better align with the latter sets of local 
data. However, a hne-tuning should be necessary for the lat¬ 
ter model to be fully applicable, in order not to overproduce 
the stringent constraints of low scatter in the loc al scaling rela¬ 
tion be tween effective radius and stellar mass (IShankar et al.l 
l2014al) . If future measurements will fully conhrm the abun¬ 
dance matching results, it may signal the need for some re¬ 
visions in key assumptions behind present hierarchical mod¬ 
els, or even the need for alternative processes beyond only 
mergers contributing to shaping the size growth o f massive 
galaxies as a function of time and environment (e.g.- lFan et al.l 
SoiOtl^ankar et al.l2013l:IStringer et al.l20l4 and references 
therein). 

The predicted sizes of centrals in groups is also systemati¬ 
cally larger by a factor of ~ 1.3 than what actually measured 
in the data, especially for the semi-empirical model character¬ 
ized by lower dynamical friction timescales. This discrepancy 


can be only partially explained by an empirical underestima¬ 
tion of the outer light profiles, as the stellar masses are well 
matched to the ones in the models, and group environments 
should be less prone to signihcant biases in s ky background 
subtra ctions, at least with respect to clusters (iBernardi et al.l 
l2013l) . From a purely theoretical point of view, this relatively 
mild discrepancy could be attributed to the progenitor disk 
sizes being too large at z ~ 1. The blue contours in the right 
panels of Fig. |4]mark the input disk sizes (Eq. [T]i for progen¬ 
itors'^ of stellar mass 0.5 and 0.7 times the remnant stellar 
mass, i.e., the typical progenitor masses in a major mergers, 
the events that in the model are responsible for creating el¬ 
lipticals (cfr. Sect. S. It is clear that the adopted disk sizes, 
although empirically motivated (Sect. S, are already notice¬ 
ably large by z ~ 1 thus clearly creating remnants that will 
inevitably overshoot the sizes of at least some ellipticals ob¬ 
served at these epochs. 

5.4. Comparison with advanced Semi-Analytic Models 

For completeness, in the same Fig . |4] we also show the 
predictions from the iGuo et all (1201 ll) semi-analytic model 
(SAM). The red, dashed lines refer to the mean properties 
of the main progenitor galaxies selected to have halo masses 
similar to those in the observational samples at z > 0 , while 
the red, dotted lines refer to the subsample of galaxies sharing 
halo masses as close as possible to those of the SDSS galaxies 
selected as described in the previous section. It is clear that 
the progenitor galaxies in the SAM are hosted by dark mat¬ 
ter halos that evolve very similarly to what predicted by the 
semi-empirical model. The stellar mass growth is instead pre¬ 
dicted to be somewhat stronger than the one in the reference 

We note that dis k sizes in Figs.CTandfsI have not been corrected by any 
projection factor (e.g.. lGonzalez et alj2009l and references therein) 































































7 




Fig. 5. — Same format as Fig.|4]but with the semi-analytic model of IGuo et all 420111) replaced by the more recent version bv IFlenrigues et al] 420141) (red 
long-dashed and dotted lines). The predicted stellar mass and size growth in the latter model is significantly lower (see text for details). 


semi-empirical model (solid lines), thus explaining why their 
mean size growth is still a factor of ~ 3 despite using /orb > 0. 
As discussed in reference to Fig.|2] models with higher values 
of /orb > 0 would be in fact expected to have less pronounced 
evolution in size at hxed progenitor masses, following from 
Eq.[3] The overall sizes at all redshifts in this SAM fall below 
the measured values. 

A more recent rendi tion of the [G uo et alj (1201 Ih model 
has been presented by iHenriaues et al.l (l2014l) . The latter 
version differs from the previous one in several key proper- 
ties. First off, they adopted the h rst-year Planck cosmology 

1 Pla nck Collaboration et al.ll2014h . at variance with lGuo et alJ 
201 ih who adopted the cosmology of the Millennium simu¬ 
lation, consistent to the one used in our semi-empirical model. 
The more relevant differences to our purposes come however 
from the baryonic sector. The num ber of mergers onto the 
centrals is signihcantly reduced in the[Henriaues et aP (1201 4) 
model with respect to iGuo et all (1201 ih . mainly because of 
two reasons. First, their flatter galaxy stellar mass function at 
the low-mass end, induces less numerous infalling satellites, 
and second, each satellite undergoes a more efficient stellar 
stripping after infall. As evident from Fig.|5] while the evolu¬ 
tion of their mean dark matter host ha los (left panels) re mains 
very similar to the one predicted by iGuo et af] (1201 Ih . their 
mean stellar mass growth is significantly reduced. Neverthe¬ 
less, centrals in clusters may still undergo sufficient stellar 
mass growth, at least when taking the SDSS group catalog as 
a local reference (open square, lower-middle panel), but defi¬ 
nitely fall short in matching the size growth (right panels). As 
expected, their me an effective radii in fac t now evolve signif¬ 
icantly less in the [Henriaues et al.l (1201 4|) mod el, at most by 
a fact or of < 1.5. Clearly setting /orb = 0 (e.g.- IShankar et all 
1201 3h and possibly reducing some of the stellar stripping in 
their model, should help increasing their predicted sizes and 
growth rates, in better agreement with the data presented in 
this work. 

6 . DISCUSSION AND CONCLUSIONS 


We have shown that by connecting high and low redshift 
BCG data via evolution of their host halo masses would imply, 
at face value, an increase since z ~ 1 of a factor ^2-3 in 
their mean stellar mass and ~ 2.5 -4 in their mean effective 
radius. Our results on the stellar mass evolution of BCGs ar e 
in line with previou s stu dies by, e.g., iLidman et alJ (1201 2l) . 
iLidman et alJ (1201 3|) . and lAscaso et al.l d2014l) . 

Although we verified that our inferred mean stellar mass 
and size growt hs hold even switch i ng to other high-re dshift 
data sets (e.g.. lAscaso et al.]l20l4 iDelave et ^l2014l) . still 
our results rely on the quality of the high redshift imaging. 
Surface brightness dimming effects, varying stellar popula¬ 
tion gradients, and/or other related spectral/photometric ef¬ 
fects, could clearly limit the general reliability of high red¬ 
shift data. Deeper imaging, possibly accompanied by accu¬ 
rate dynamical measurements, will certainly be able to pro¬ 
vide more stringent constraints on the actual stellar mass and 
size growths of BCGs. 

At face value, we And a good match between 
observationally-derived stellar mass and size growth 
rates, and those predicted from state-of-the-art, hierarchical 
galaxy evolution models. In the specific, our best models 
predict an average number of ^ 0.5 -0.9 major mergers and 
^7-8 minor mergers since z ^ 1 (Fig. B, with mass ratios 
fi in the progenitors of 0.3 < p < 1.0 and 0.01 < /r < 0.3, re- 
spectiv ely. In accordance to our re sults, in the ALFIAMBRA 
survey lUopez-Sanjuan et al] (l2014l) recently estimated in the 
same redshift interval an average number of major mergers 
per gala xy of Nm = 0-57 ib 0.05 fo r red luminous galaxies 
(see also iLopez-Saniuan et al.ll2012l) . Our predicted ayerage 
major merger ra tes of R ^ (0.07-0.1 l)Gyr~* is also in good 
agreement with iLackner et alJ (l2014h inferred from a mass 
completed sample of massiye galaxies R ^ 0.1Gyr“/ down 
to R ^ 0.04Gyr~* for quiescent galaxies alone. 

iLiu et alJ (1201 5h recently inferred from BCG deep, high- 
resolution imaging a major merger rate of (0.55 ± 0.27) Gyr“' 
at z 0.43, concluding that present-day BCGs should haye 
increased their stellar mass by ^ (35 ± 15)% yia major dry 
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mergers since z = 0.6, in good agreement wi th what de¬ 
rived from our best models. Ilnagaki et al.l (1201 5h . in analogy 
to our methodology, constructed from the Planck Sunyaev- 
Zel’dovich cluster catalog galaxy progenitor-descendant pairs 
of BCGs, finding an average 2-14% growth in the redshift 
range 0.2-0.4 (see also lOliva-Altamirano et al.ll2014h . in good 
agreement with our average growth of ~ 10% in the same 
redshift interval. Clustering measureme nts provide aver¬ 
age merger rates of ~ 0.024Gyr“' (e.g., IWake et ^ 120081 : 
IWhite et al.l2008l and references therein) at z < 0.6 somewhat 
lo wer than our estim ates. 

iMan et aP (1201 4l) from a mass-complete sample of massive 
galaxies from UltraVISTA/COSMOS, complemented with 
deeper, higher resolution SDHSTh-CANDELS data, have also 
estimated a major merger rate R ~ 0.1 Gyr“' for comparably 
massive galaxies, at least since z 1.5 (e.g., their Fig. 12). 
The authors suggest that such an amount of major merging 
alone could be sufficient to explain the observed number den¬ 
sity evolution for massive galaxies (logMstar/440 > 11.1). In¬ 
deed, although still debated, the combined information from 
current stellar mass, clustering, and weak lensing measure¬ 
ments on the stellar mass and environments of intermediate- 
redshift BCGs tend to favor a noticeable increase in the num¬ 
ber d ensity of the most mas sive galaxies (e.g., Buitrago et al. 
2013. l^ankar et aT]l2014bl an d references therein). 

However. iMan et al.l (1201 4l) also measure a minor merger 
rate of R ^ 0.07 - 0.1 Gyr“* that would be insufficient by 
itself to explain the rapid mean size evolution of a fac¬ 
tor of ~ 2 -2.5 since z 1 as inferred f r om our data (see 
also, e.g., iHuertas-Comnanv et al.l 1201 3bt iMan et akl 1201 4l 
van der Wei et al.ll2014il . However, we should note that the 
Man et al.l (1201 4|) estimate is limited to minor mergers with a 
mass ratio 0.1 < /r < 0.3, while hierarchical models predict, as 
in the analysis presented in this work, that extending this limit 
to 0.01 < /r < 0.3, i ncreases the rate by a factor of > 5 (e.g., 
iHopkins et al.ll20f^ and references therein). Such a signifi¬ 
cantly increased rate of minor mergers has a relatively lower 
repercussion on the increase of the stellar mass , but more ef¬ 
fectively impacts the o verall size growth (e.g.. iTruiillo et all 
l2007HNaah et~^l2009h . 

We showed that the models with lower dynamical friction 
timescales, though disfavored on more general grounds be¬ 
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caus e predicting too loose scaling relations a t low redshifts 
(e.g., iBower et^l2006l IShankar et al.ll2014^ . may still be a 
viable solution for properly growing massive, central galaxies 
in cluster-sized environments. The number of predicted ma¬ 
jor mergers in the latter model is also similar to the one pre¬ 
dicted by the reference model with longer dynamical friction 
timescales, though the number of minor mergers is double. 
Unfortunately, present data cannot still reliably constrain the 
number of minor mergers to distinguish among these models. 

To conclude, in order to limit the effect of progenitor bias, 
we have gathered a number of data on high redshift bright¬ 
est group and cluster galaxies, evolved their mean host halo 
masses down to z = 0 along their merger trees main progenitor 
branches, and assigned as their descendants local SDSS cen¬ 
tral galaxies matched in host halo masses. At face value, the 
direct comparison of local “descendants” and z ^ 1 “progeni¬ 
tors”, suggests an increase in mean stellar mass of a factor of 
> 2, and a relative increase of mean effective radius of > 2.5. 
We found that the most favored hierarchical models, with no 
strong stellar stripping and null orbital energies, yield a factor 
of ~ 2.5 - 3 growth in median size, in good agreement with 
what inferi'ed from the matching between high and low red¬ 
shift data, especially for the stellar mass growth of centrals in 
group-sized, high-redshift halos. 
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